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1)  Spectral  methods  for  discontinuities. 

We  have  showed  that  for  linear  problems  spectral  methods  yield  highly 
accurate  information  even  for  problems  with  discontinuities.  Appendix  A 
is  a  summary  of  invited  lectures  given  by  Dr.  Gottlieb  at  the  Conference 
on  Numerical  Methods  in  Fluid  Dynamics  held  at  Reading,  England  in  March 
1985. 

Basically,  we  showed  that  for  linear  problems  the  Fourier  or  Fourier- 
Chebyshev  coefficient  of  the  solution  are  obtained  with  spectral 
accuracy.  We  also  showed  several  techniques  to  recover  spectrally 
accurate  pointwise  values  from  the  Fourier  coefficients  of  discontinuous 
data.  We  have  simulated  also  the  wedge  flow  using  the  full  Euler 
equations  in  2D  and  found  that  also  in  this  case  the  results  contain 
spectrally  accurate  information. 
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We  have  made  progress  in  three  areas 


2)  Boundary  conditions  for  incompressible  flows. 

We  have  developed  simple  techniques  to  achieve  high  order  accurate  time 
integrations  for  incompressible  flow  equations.  The  results  are 
summarized  in  Apper  !ix  B. 
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3)  We  have  carried  out  benchmark,  testing  of  large  hydrodynamics  code  on 
supercomputers . 


Spectral  codes  have  been  benchmarked  on  a  large  variety  of  American  and 
Japanese  supercomputers,  including 


Cray  Research  - 
Control  Data  - 
Hitachi  - 
Fuj i tsu  - 
NEC  - 


Cray  IS,  Cray  XMP24,  Cray  2 
Cyber  205-642 
Hitachi  SA810/20 
Fujitsu  VP-200,  VP-400 
NEC  SX-2. 


We  have  found  that  spectral  hydrodynamics  codes  are  highly  vector izable  • 
They  rely  heavily  on  mu  1 1 i-d iraens ional  fast  Fourier  transforms  (FFTs). 

The  following  performance  was  recorded: 

SHEA£3  -  Free-slip  convection/shear  flow  code  -  64^  resolution 


Cray  IS 

80 

MFlops 

(i 

processor) 

Cray  XMP24 

260 

MFlops 

(2 

processors ) 

Cray  2 

1100 

MFlops 

(4 

processors ) 

Cyber  205 

210 

MFlops 

(1 

processor) 

Fujitsu  VP-2D'J 

140 

MFlops 

(1 

processor) 

Fujitsu  VP-4i; 

180 

MFlops 

(1 

processor) 

Hitachi  SA810/J 

110 

MFlops 

(1 

processor ) 

NEC  SX-2 


380  MFlops  (1  processor) 
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Novel  programming  is  sometimes  necessary  to  take  maximal  advantage  of  some 
of  the  newer  architecture  machines.  For  example,  on  the  Cray-2,  it  is 
necessary  to  make  efficient  use  of  the  local  16K.  memories  on  each 
processor  to  achieve  two  simulataneous  paths  to  memory  and,  hence,  achieve 
full  speed  on  FFTs. 


